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We present a theoretical analysis of the appearance of entanglement in non-interacting mesoscopic 
structures. Our setup involves two oppositely polarized sources injecting electrons of opposite spin 
into the two incoming leads. The mixing of these polarized streams in an ideal four-channel beam 
splitter produces two outgoing streams with particular tunable correlations. A Bell inequality test 
involving cross-correlated spin-currents in opposite leads signals the presence of spin-entanglement 
between particles propagating in different leads. We identify the role of fermionic statistics and 
projective measurement in the generation of these spin-entangled electrons. 



Quantum entangled charged quasi-particles are per- 
ceived as a valuable resource for a future solid state based 
quantum information technology. Recently, specific de- 
signs for mesoscopic structures have been proposed which 
generate spatially separated streams of entangled parti- 
cles 0,0, SB- In addition. Bell inequality type measure- 
ments have been conceived which test for the presence 
of these non-classical and non-local correlations 0, 0|. 
Usually, entangled electron-pairs are generated through 
specific interactions (e.g., through the attractive inter- 
action in a superconductor or the repulsive interaction 
in a quantum dot) and particular measures are taken to 
separate the constituents in space (e.g., involving beam 
splitters and appropriate filters). However, recently it 
has been predicted that non-local entanglement as sig- 
nalled through a violation of Bell inequality tests can be 
observed in non- interacting systems as well 5, (I 0,13 ■ 
The important task then is to identify the origin of the 
entanglement; candidates arc the fermionic statistics, the 
beam splitter, or the projection in the Bell measurement 
itself 0,0]. 

Here, we report on our study of entanglement in a non- 
interacting system, where we make sure, that the parti- 
cles encounter the Bell setup in a non-entangled state. 
Nevertheless, we find the Bell-inequality to be violated 
and conclude that the concomittant entanglement is pro- 
duced in a wave function projection during the Bell mea- 
surement. We note that wave function projection as a 
resource of non-local entanglement is known for single- 
particle sources (Fock states) 0|, a scheme working for 
both bosons and fermions. What is different in Refs. 
0000 and in the present work is that the sources are 
many-particle states in local thermal equilibrium. It is 
then essential that one deals with fermions; wave function 
projection cannot create entanglement out of a thermal 
source of bosons. 

The generic setup for the production of spatially sep- 
arated entangled degrees of freedom usually involves a 
source injecting the particles carrying the internal degree 
of freedom (the spin 0000] or an orbital quantum 
number 0|) and a beam splitter separating these 
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FIG. 1: Mesoscopic normal-metal structure with a beam split- 
ter generating two streams of electrons with tunable correla- 
tions in the two outgoing arms 'u' and 'd'. The source (left) 
injects polarized (along the z-axis) electrons into the source 
leads 's' and 's'. The beam splitter mixes the two incoming 
streams with a mixing angle i9. The scattered (or outgoing) 
beams are analyzed in a Bell type coincidence measurement 
involving spin-currents projected onto the directions ±a (in 
the 'u' lead) and ±b (in the 'd' lead). The injection reser- 
voirs are voltage (V) biased against the outgoing reservoirs. 
The Bell inequality test signals the presence of entanglement 
within the interval |i?-45°| < 12.235°. We relate this entan- 
glement to the presence of spin-triplet correlations in the pro- 
jected part of the scattered wavefunction describing electron- 
pairs distributed between the arms. 

particles in space, see Fig.^ In addition, 'filters' may be 
used to inhibit the propagation of unwanted components 
into the spatially separated leads [Jl , thus en- 
forcing a pure flow of entangled particles in the outgoing 
leads. The successful generation of entan glem ent then is 
measured in a Bell inequality type setup [lOj . A surpris- 
ing new feature has been recently predicted with a Bell 
inequality test exhibiting violation in a non-interacting 
system 0000; the question arises as to what pro- 
duces the entanglement manifested in the Bell inequality 
violation and it is this question which we wish to address 
in the present work. In order to do so, we describe the- 
oretically an experiment where we make sure, that the 
particles are not entangled up to the point where the 
correlations are measured in the Bell inequality setup; 
nevertheless, we find them violated. We trace this vi- 
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olation back to an entanglement which has its origin in 
the confluence of various elements: i) the Fermi statistics 
provides a noiseless stream of incoming electrons, ii) the 
beam splitter mixes the indistinguishable particles at one 
point in space removing the information about their ori- 
gin, in) the splitter directs the mixed product state into 
the two leads thus organizing their spatial separation, 
iv) a coincidence measurement projects the mixed prod- 
uct state onto its (spin-)entangled component describing 
the electron pair split between the two leads, v) mea- 
suring the spin-entangled state in a Bell inequality test 
exhibits violation (the steps iv) and v) are united in our 
setup). Note, that the simple fermionic reservoir defin- 
ing the source in Ref.0 injects spin-entangled pairs from 
the beginning; hence an analysis of this system cannot 
provide a definitive answer on the minimal setup pro- 
viding spatially separated entangled pairs since both the 
source and/or the projective Bell measurement could be 
responsible for the violation. 

Below, we pursue the following strategy: We first de- 
fine a particle source and investigate its characteristic 
via an analysis of the associated two-particle density ma- 
trix. We then define the corresponding pair wave func- 
tion (thus reducing the many body problem to a two- 
particle problem) and determine its concurrence follow- 
ing the definition of Schliemann et al. for indistin- 
guishable particles (more generally, one could calculate 
the Slater rank of the wave function, cf. Ref.^3 here, we 
deal with a four-dimensional one-particle Hilbert space 
where the concurrence provides a simple and quantita- 
tive measure for the degree of entanglement). For our 
specially designed source we find a zero concurrence and 
hence our incoming beam is not entangled. We then go 
over to the scattering state behind the (tunable) beam 
splitter and reanalyze the state with the help of the two- 
particle density matrix. We determine the associated 
two-particle wave function and find its concurrence; com- 
paring the results for the incoming and scattered wave 
function, we will see that the concurrence is unchanged, 
a simple consequence of the unitary action of the beam 
splitter. However, the mixer removes the information 
on the origin of the particles, thus preparing an entan- 
gled wave function component in the output channel. 
Third, we analyze the component of the wave function 
to which the Bell setup is sensitive and determine its 
degree of entanglement; depending on the mixing angle 
of the beam splitter, we find concurrencies between 
(no entanglement) and unity (maximal entanglement). 
Finally, we determine the violation of the Bell inequal- 
ity as measured through time-resolved spin-current cross- 
correlators and find agreement between the degree of vi- 
olation and the degree of entanglement of the projected 
state as expressed through the concurrence. 

Our source draws particles from two spin-polarized 
reservoirs with opposite polarization directed along the 
z-axis. The polarized electrons are injected into source 



leads 's' and 's' and are subsequently mixed in a tun- 
able four-channel beam splitter, see Fig. ^ The outgoing 
channels are denoted by 'u' (for the upper lead) and 'd' 
(the 'down' lead). The spin-correlations in the scatter- 
ing channels 'u' and 'd' are then analyzed in a Bell in- 
equality test. The polarized reservoirs are voltage biased 
with eV = i1bH/2 equal to the magnetic energy in the 
polarizing field H; the incoming electron streams then 
are fully polarized (the magnetic field is confined to the 
reservoirs). 

The spin-correlations between electrons in leads 'x' and 
'y' are conveniently analyzed with the help of the two- 
particle density matrix (or pair correlation function) 

gf{x, y) = Tt{p^I, {x)^,^ {y)^y,, (y)^,,, (x)) (1) 

with trace over states of the Fermi sea. Here, are 
field operators describing electrons with spin a in lead 
'x' and p is the density operator. The pair correlation 
function is conveniently expressed through the one- 
particle correlators G^^^{x,y) = ('I'^^(a;)5'yg-(y)), 

97 i^^y)- /o)G7^3(o) - G7a^~y)Gi-..Sy~^) ■ (2) 

The one-particle correlators can be written in terms 
of a product of orbital- and spin parts, G7six,y) = 
G^y{x,y)x^^{'^,<^), and spht into equihbrium and excess 
terms, 

Gll{x,y) = Gcq{x,y)x^^{(J,a) + Gcx{x,y)Xcl{<^,^), 

(3) 

with Gcx{x,y) vanishing at zero voltage V and zero po- 
larization field H. 

In order to find the two-particle density matrix in the 
source leads 's', 's' we make use of the scattering states 

^-, = ^ e''="afe<, + e-^'=^(cosz?e-''^Cfe, + siniJe^'^dfc,), 

ka 

^, = ^e''=-6,, + e-^'=^(cos^9e^'^dfc„ - sini9e-^'^Cfc„), 

kcr 

where aka , bka denote the annihilation operators for elec- 
trons in the source reservoirs s and s with momentum k 
and spin cr Gti i polarized along the z-axis and time evo- 
lution cx exp{—iekt/h), = h^k'^/2m; the operators Cka 
and dka annihilate electrons in the reservoirs attached 
to the outgoing leads 'u' and 'd', respectively. Also, we 
make use here of the standard parametrization of a re- 
flectionless four-beam splitter, 

f n\ _ f e'^^cos^ -e^'^sinz? \ / s \ 

l^dy l^g-^'^sini? e-*^cosi9 ) \ s ) ' ^ ' 

with the angles i) G (0,7r/2), (p,tp £ (0, 27r); without 
loss of generality we will assume ip = -ip — in what 
follows. The orbital part of the one-particle correlator 
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G^^{x — y) = G{x — y) takes the form 

Gcq(a;) ^ , (5) 

nx 

with kv ~ k-p{eV/ep) and (kp) the Fermi energy (wave 
vector) in the unbiased system. The spin factors for the 
equiUbrium and excess parts read, 

= (7) 
Xllia, a) = {a\ T) (T |a) , ^{a, a) ^ (a| [) (i 

the latter describing the injection of polarized electrons 
into the leads 's' and 's'. Finally, the cross correlation 
function between the source leads vanishes, G^%{x — y) = 
0, and the final result for the excess part of the pair 
correlation function between source leads reads 

W{x,y)]^^ = \GU^)\^a,\]){]\a,){a^\i){i\a,). (8) 

This result then describes the injection of two uncor- 
related streams of polarized electrons into the leads 's' 
and 's'. Furthermore, statistical analysis |0| tells that 
the Fermi statistics enforces injection into each lead of 
a regular stream of particles separated by the single- 
particle correlation time ry = Ti/eV. The full many 
body description then is conveniently reduced to a two- 
particle problem where the two reservoirs inject a se- 
quence of electron pairs residing in the wave function 
■^11 = - 0i^0^^]/V2 with </.sT iM) the single- 

particle wave functions associated with electrons in the 
upper (lower) source lead. This wave function is a simple 
Slater determinant and hence non-entangled according 

to [13]. 

Next, we extend the above analysis to the outgoing 
leads 'u' and 'd'. The scattering states in the outgoing 
leads take the form 

^'u =^e"'''"^Cfe<T+e*''="(cosi^afc<, - sini?6fc<,), 

k(7 
ka 

The excess particles injected by the source leads now are 
mixed in the beam-splitter and thus non- vanishing cross 
correlations are expected to show up in the leads 'u' and 
'd'. The one-particle correlation function assumes the 
form Q with the orbital correlators |(SJ| and ^ and spin 
correlators 

X^^(a,a) = (a|a), x e V, 'd', (9) 
X:^^ia, a) = cos2 t) (T \a) + sin^ i) (i \a) , 
Xt^{<y, a) = sin^ ]) {] \a) + cos^ §{a\ [) {[ \a) , 

= cosi?sin^[(a|T)(T|a)-(aU)a|a)]. 



Evaluating the excess part of the two-particle cross- 
correlations between the leads 'u' and 'd' at the sym- 
metric position X = y we find 

[gfix,x)]^^=\G,^iOr (10) 
[cos4z9(ai|t>(Tk4)(rT2U)(i|a3) 

-Ksin4i?((TiU)(i|a4)(a2|t)(Tk3) 

+ cos2 ^ sin^ ^(ai | T) (T ks) {<J2\l) (i 1^4) 

+ cos2 ^ sin2 ^(ai | i) (i jag) {'^2 \ T) (T ^4)] ■ 

Hence, a symmetric splitter (i? — 7r/4) produces the 
spin correlations of a triplet state [\xtf) = I T)u| i)d+ 
I i)u| T)d]/'\/2 involving two electrons separated in differ- 
ent leads 'u' and 'd' but at equivalent locations x — y. 
The general case with arbitrary mixing angle •& results in 
a density matrix describing a pure state involving the su- 
perposition \xtf) + cos2i?|xsg ) of the above triplet state 
and the singlet state [\x'^^) = | t)u| i)d - | i)u| T)d]/\/2. 
The analoguous calculation for the two-particle density 
matrix describing electrons in the same outgoing lead 'x' 
equal 'u' or 'd' points to the presence of singlet correla- 
tions, 

[grix,y)],^ - \GUO)\^<jM{<j,\as) (11) 
- \Gcxix - y)f((Ji\a:i)(a2\(T4). 

Again, the above results can be used to reduce the 
problem from its many-body form to a two-particle prob- 
lem. Given the incoming Slater determinant we 
obtain the scattered state VP^^t through the transfor- 
mation (psi — > cosi}(/)u^ + sini?(/)d| describing scattered 
spin-t electrons originating from the source lead 's' and 
0si — > — sini9 0ui +coai)(f)di for excess spin-| electrons 
from 's' (the wave functions (p^a — i'xXa describe elec- 
trons with orbital (spin) wave function 0x (Xo-) propagat- 
ing in lead 'x'). The resulting scattering wave function 
has the form 

= sin^cos^[0:,Cx^^ - <l^'^4>lxll] 
+<I>-dXt"-)-cos2^$-dXsg, (12) 

where the first two terms describe the propagation of a 
spin-singlet pair with the wave function — {x\x\ ~ 

XjXp/v^ in the upper and lower lead. The last two 
terms describe the component where the electron pair is 
split between the 'u' and 'd' leads; it is a superposition 
of singlet- and triplet states [xtl = (X|Xj + x\x^])/"^) 
with corresponding symmetrized and anti-symmetrized 
orbital wave functions i>Jjj = ((/)u(?!>^ + (/'Ji(?!'u)/2 and = 
{'t'n^'d ~ 'Pd'f'u)/'^- The entanglement present in these 
wave functions is easily determined using the formalism 
developed by Schliemann et al. |l ij : The wave func- 
tion associated with a pair of electrons can be written in 
terms of a single-electron basis {(f>i}, '^^^ = ^l^ij^''] 
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where the anti-symmetric matrix Wij = —Wji guaran- 
tees for the proper symmetrization. The analysis sim- 
phfies drastically for the case where the one-particle 
Hilbert space is four-dimensional; then the concurrence 
C{'^) = Sy'detwi^ gives a quantitative measure for the 
entanglement present in the wave function C(5') = 
for a non-entangled state and C{^) = 1 for a fully entan- 
gled wave function. For our setup the one-particle basis 
is defined as {(/)u]', '/'ui, 0dT, (f'di} and the matrix ?i'(5'out) 
describing the scattered state H12|l assumes the form 



via the current difference correlators 



1 

V2 





sin 2'0/2 


— COS^ 1? 



-sin2i?/2 


- sin^ 




cos^i? 

sin^ 

sin2?9/2 

-sin2i?/2 



The concurrence of the scattering state ()12f) vanishes, 
hence v[/°ut is non-entangled and takes the form of an el- 
ementary Slater determinant. Next, let us analyze the 
concurrence of that part of the scattering wave func- 
tion to which our coincidence measurement in leads 'u' 
and 'd' is sensitive. The component describing the two 
particles split between the leads reads 5*^^^ — <&udXtr + 
cos2i?3>Jj^^Xsg' '^f- 02l- This projected state is described 
by the matrix 



1 

71 






- sin^ 
cos^ i9 





sin^i? 





COS^ 1? 







from which one easily derives the concurrence — 
sin^ 2??; we conclude that the component detected in 
a coincidence measurement is entangled. Furthermore, 
the concurrence is equal to unity for the symmetric split- 
ter = 7r/4 where we deal with a maximally entangled 
triplet state (note the loss of information about which 
electron (from 's' or 's') enters the lead 'u' or 'd'). We 
conclude that a Bell inequality test sensitive to the split 
part of the wave function will exhibit violation. We at- 
tribute this violation to the combined action of the 
splitter where the information on the identity of the par- 
ticles is destroyed and the entangled component '^^^^ is 
'prepared' and ii) the wave function projection inherent 
in the coincidence measurement and 'realizing' the en- 
tanglement. 

The Bell type setup ^|] in Fig. ^ measures the cor- 
relations in the spin-entangled scattered wave function 
^'q^^. It involves the finite-time current cross-correlators 
Ca,bix,y;T) = ((/a(a:,T)/b(y,0))) between the spin- 
currents Ia{x,T) projected onto directions a (in lead 
'u') and partners lb{y,0) (in lead 'd') projected onto b. 
These correlators enter the Bell inequality (a and b de- 
note a second set of directions) 

\E{a,h) - E{a.,h) + E{a,h) + E{a,h)\ < 2 (13) 



E{a,h) 



([/a(r)-/-a(r)][/b(0)-/-b(0)]) 

([/a(r)+/-a(r)][/b(0)+/-b(0)])' 



(14) 



The cross-measurement in different leads implies that 
the setup is sensitive only to the spin-entangled split- 
pair part of the scattering wave function and hence 
the Bell inequality can be violated. Making use of the 
field operators and \l/d describing the scattering states 
in the outgoing leads, we determine the irreducible cur- 
rent cross correlator and factorize into orbital and spin 
parts, Ca,b(a;,2/;T) = Ca;,y(T)i^a,b, with Fa,b accounting 
for the spin projections. Using standard scattering the- 
ory of noise [iSj, one obtains the orbital cross-correlator 
(only the excess part gives a finite contribution) 



Cx,y{T) 



2 



sin^2i? 



eV{T- 



■a{T-T-,9), (15) 



with a{T,0) = TT'^9^/smh^[Tr9T/n], t± = {x±y)/v^, 9 
the temperature of the electronic reservoirs, and Wp the 
Fermi velocity. In order to arrive at the result (|15|) we 
have dropped terms small in the parameter \e' — e\/e-p 
[l^ . The spin projection Fa,b assumes the form 

Fa,b = (a| T) (T |b) (b| T) (T |a) + (a| j) (i |b) (bji) (j |a) 
- (a| t) (T |b) (bli) a |a) - (ali) (i |b) (b| T) (T |a) . 

We express this result in terms of the angles 9^ and ip^ 
describing the direction of magnetization in the 'u' lead 
filters and 9b, (pb referring to the filters in the 'd' lead and 

find that Fa,b = i^-a,-b = Fatb- ^-a,b = i^a,-b = i^a^b 

and 

F^y^ — {1± cos 0a cos ^bTcos ifab siu 0a sin 0b)/2, 

with ipab = fa ^ fb- The correlator ^^(a, b) takes the 
form 



E{ai,h) 



^C^Ar)[F+b-Pa,b]+^- 



with A± = [(/a) ± (/-a>][(/b) ± (I-b)]- Evaluat- 
ing the projected current averages one obtains A_ = 
-e2(2el///i)2cos6'aCOse'bCOs2 2d and A+ = e^{2eV/h)^. 
The triplet state is rotationally invariant within the plane 
^a = ^b = 7''/2 and choosing filters within this equatorial 
plane the Bell inequality takes the form 



Cx,y (r) [cos l/?ab - cos yJab + cos yab + COS (/?ab] 



2Cx^yiT)+A^ 



< 1. 



Its maximum violation is obtained for the set of angles 
'/'a 0, (/3b = 7r/4, (/^a = 7r/2, ipb = 37r/4, 



Ebi — 



2CxJt)+A+ 



< 



1 



(16) 
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Evaluating the above expression in the limit of low tem- 
peratures 9 < eV and at the symmetric position x = y, 
we arrive at the simple form 

sin2 2i9 sin^feyW^i) 1 , ^ 
!^ —L < — (17) 

2{eVT/n)'^ -sm'^2^sm^{eVT/h) - V2' ^ ' 

We observe that the violation of the Bell inequality is 
restricted to short times r < t^i = Ty = h/eV 
the relevance of a coincidence measurement involving the 
short time ry was noticed in Refs. 0,0]). Furthermore, 
the violation strongly depends on the mixing angle d of 
the beam splitter, with a maximal violation realized for 
a symmetric splitter d = n/A generating a pure triplet 
state across the two arms. The Bell inequality cannot be 
violated for asymmetric splitters with \d~ 7r/4| > 0.2135 
(corresponding to an angular width — 45°| > 12.235°): 
evaluating the BI 1)17(1 at zero time difference (i.e., in a 
coincidence measurement) we find the condition 

sin^2i? 1 
9 < 18) 

from which one derives the critical angle = 
[arcsin(2/(v^+l))i/2]/2 = 0.572 (or = 32.765°). The 
appearance of a critical angle naturally follows from the 
fact that the measured wave function component 
assumes the form of a simple Slater determinant in the 
limits d = 0, 7r/2 and hence is not entangled. Note that 
the product of average currents A+ is the largest term in 
the denominator of (|16|l and hence always relevant. 

In conclusion, we have described a mesoscopic setup 
with a source injecting non-entangled electron pairs into 
two source leads 's' and 's'. Subsequent mixing of these 
particle streams in a four-channel beam splitter does not 
generate entanglement between the particles in the two 
output leads 'u' and 'd'. However, proper mixing of the 
incoming beams in the splitter removes the information 
on the path of the incoming particles and generates a 
wave function component describing electrons split be- 
tween the leads 'u' and 'd' which is entangled. It is this 
component which manifests itself in the coincidence mea- 
surement of a Bell inequality test and proper violation is 
observed at short times. This analysis answers the ques- 
tion regarding the origin of entanglement observed in the 
Bell inequality test applied to the present non-interacting 
system. A modified setup where the particles propagate 
downstream after a coincidence measurement lends itself 
as a source for spin-entangled particles, cf. Ref. |3 

It is interesting to analyze the setup described in Ref.Ja 
in the light of the findings reported here. The setup in p| 
involves a simple normal reservoir injecting pairs of elec- 
trons into a source lead which are subsequently separated 
in space by a beam splitter. The injected pairs reside in 
a spin-singlet state involving the identical orbital wave 
function, '^['^ = (j)l(j)lxlll the entanglement observed in a 
Bell inequality test then has been attributed to the en- 
tanglement associated with this spin-singlet state. One 



may criticise, that this incoming singlet, being a simple 
Slater determinant, is not entangled according to the def- 
inition given by Schliemann et al. However, after 
the beam splitter the orbital wave function (f>s is delo- 
calized between the two leads, (jig ^ ^ — isu^u + isd^d, 
with tsn and tgd the corresponding scattering amplitudes. 
While the scattered state remains a Slater determinant 
^out = ^^^^xll, the singlet correlations now can be ob- 
served in a coincidence measurement testing the cross- 
correlations between the leads 'u' and 'd'. Hence the 
spin-entanglement is produced by the reservoir, but its 
observation requires proper projection. It is then difficult 
to trace a unique origin for the entanglement manifested 
in the violation of a Bell inequality test. The appropriate 
setup to address this question should involve a reservoir 
injecting particles with opposite spin residing in a Slater 
determinant of the form 'I'J^ = — (j)l^(j)l.^]/^/2 
which is not entangled in the spin variable. Such an 
analysis has been presented here with the result, that 
the orbital projection in the coincidence measurement is 
sufficient to produce a spin-entangled state. 
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